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Abstract 

In a previous work, assuming that the nucleus can be treated as a perfect fluid, we have stud- 
ied the propagation of perturbations in the baryon density. For a given equation of state we 
have derived a Korteweg - de Vries (KdV) equation from Euler and continuity equations in non- 
relativistic hydrodynamics. Here, using a more general equation of state, we extend our formalism 
to relativistic hydrodynamics. 
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I. INTRODUCTION 

Long ago [l| it was suggested that Korteweg - de Vries solitons might be formed in the 
nuclear medium. In a previous work [2] we have updated the early works on the subject 
introducing a realistic equation of state (EOS) for nuclear matter. We have found that these 
solitary waves can indeed exist in the nuclear medium, provided that derivative couplings 
between the nucleon and the vector field are included. These couplings lead to an energy 
density which depends on the Laplacian of the baryon density. For this class of equations of 
state, which is quite general (as pointed out in j^, 4|), perturbations on the nuclear density 
can propagate pulse without dissipation. 

During the analysis of several realistic nuclear equations of state, we realized that, very 
often the speed of sound Cg is in the range 0.15 — 0.25. Compared to the speed of light these 
values are not large but not very small either. This suggests that, even for slowly moving 
nuclear matter, relativistic effects might be sizeable. This concern motivates the extension 
of the f„n„*. p..ese,.ed i„ Q .e.a«.ic M.ody„^^^^^ 

In the next section we review the most relevant equations writting them in an appropriate 
form for the subsequent manipulations. In the following section we discuss three models for 
the equation of state and in the next section we derive the KdV equation for the proposed 
models and present their solutions. In the final section we present our conclusions. 

II. RELATIVISTIC HYDROYNAMICS 

In this section we review the main expressions of relativistic hydrodynamics. In natural 
units (c = 1) the velocity four vector u'^ is defined as: 

= {u',u) ={i,^v) (1) 

where 7 is the Lorentz contraction factor given by: 

i = {i-vr'^' (2) 

The velocity field of matter is v = v{t,x,y,z) and thus u'^Uiy = 1. The energy-momentum 
tensor is, as usual, given by: 

T^Lu = {£ + p)u^,Uy - pg^u (3) 
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where e and p are the energy density and pressure respectively. Energy-momentum conser- 
vation is ensured by: 

= (4) 

The projection of jU]|_pnto a direction perpendicular to gives us the relativistic version 
of Euler equation 



dv -> 1 / -> dv\ 

The relativistic version of the continuity equation for the baryon number is 

O^jb" = (6) 
Since Jb" = u'^Pb the above equation reads 

^^{PB7)+V-{PB1V)=0 (7) 



The enthalpy per nucleon is given by [6[: 

dh = Tds + Vdp (8) 

where V = 1/pb is the specific volume. For a perfect fluid {ds = 0) the equation above 
becomes dp = pBdh and consequently: 

^ ^, 9p dh . . 

Vp = PbVA. - = p.^ (9) 

Inserting ([9]) in ([5]) we find: 

dv ,^ -L^ ^ pb -^dh\ . 



Recalling the Gibbs relation 



e-Ts + p = pBpB (11) 
and considering the case where T = we obtain: 

e + p = PbPb (12) 

where e, p, pb and are the energy density, pressure, baryochemical potential and baryon 
density respectively. Inserting ([21) and ( fT2l) into ( ITOi) we obtain: 

_ + („.V)„=-L— i(vA + .^) (13) 
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We close this section comparing the relativistic and non-relativistic versions of the Euler 
and continuity equations. The latter were presented in 

^^"^ + V ■ (PBV) = (14) 



^+{v-V)v = -[-Wh (15) 



dt 

dt ^^'^ ''^ \M 
and the former are ([7]) and f|T3l) : 

^(Pi?7) + V-(PB7^) = (16) 

_ + („.V).= -^(v/> + »-) (17) 

The two pairs are similar. The differences are only in the 7 factors and in the last term 
of (fT7j) . where the appearance of time derivative reflects the symmetry between space and 
time. Since the enthalpy per nucleon may also be written as [2|, 

de 
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(18) 



dpB 

it becomes clear that the "force" on the right hand side of (flSj) and ( |T7I) will be ultimately 
determined by the equation of state, i.e., by the function s^Pb)- 



III. EQUATION OF STATE 

Equations f|T5l) and f|T7|) contain the gradient of the derivative of the energy density. If e 
contains a Laplacian of Pb-, i.e., e oc ... + ...V'^Pb + then fllSI) and (11 7p will have a cubic 
derivative with respect to the space coordinate, which will give rise to the Korteweg-de Vries 
equation for the baryon density. The most popular relativistic mean field models do not 
have higher derivative terms and, even if they have at the start, these terms are usually 
neglected during the calculations. 

f] rn 

In |2l] we have added a new derivative term to the usual non- linear QHD ^] , given by 

Cm = -^^id^d^V,)^^ (19) 

where, as usual, the degrees of freedom are the baryon field ip, the neutral scalar meson field 
and the neutral vector meson field V^, with the respective couplings and masses. The new 
term is designed to be small in comparison with the main baryon - vector meson interaction 
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term g^ip'-j^V^'il). Folowing the standard steps of the mean field formahsm we have arrived 
at the following expression for the energy density 



9v 



pI + 



(M* - M) 



9s 

2 



+ 



V 



(2vr) 



kp _, h 

^9s 



(20) 
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where rj is the baryon spin-isospin degeneracy factor, M* stands for the nucleon effective 
mass (given by M* = M — gs<f)o) and the constants b, c, Qs and were taken from 10|. 
Although Eq. (l20i) was obtained with the help of a specific Lagrangian taken from 10| and 
a prototype Laplacian interaction f|T9l) . the above form of the energy density follows quite 
naturally from an approach based on the density functional theory ^ , for a wide variety of 
underlying Lagrangians. 

We now follow the treatment developed in to obtain the Korteweg-de Vries 

equation in one dimension through the combination of ( fT6l) and ( |T71) . With the help of ( 120|1 
we first calculate the energy per nucleon given by E = e/ps- We next perform a Taylor 
expansion of E around the equilibrium density po up to second order: 



Eipn) = Eipo) + - 



1/d^E 



2 \dpB' 



(pB - Pof 



PB=PO 



where the first order term vanishes because of the saturation condition: 



d 



-M 



dpB \PB 

We arrive at (for more details see ^ 
E{pb) = E{po) + 



PB=PO 



9v 



/dE_ 



rrir 



B / PB=PO 



IMc 



{"^ Pb) + -^{PB - PoY 



Po 



The enthalpy per nucleon may also be written as 9|: 

dE 



h = E + pB- 



Using 



and 



to evaluate 



dh 
dt 



dpB 

and its derivatives we find: 
2Mc 2 



:^PBypB 

Po^ 



9v 



Po 



^PB + ^V(V>b) 



my-" 



SMc,"^ dpB 2Mcs'^dpB , c/^^ d 



Po^ ot 



Po 



dt 



+ 



m^*^ dt 



(V'pb) 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



The expressions ( 125|) and (126|) . from now on referred to as model I, will be inserted into 
f|T7|) as it will be seen in the next section. 

We shall now consider a more general expression for the energy density 

e = aipB + «2Pb^ + «3Pb^ + Ppb'^^Pb (27) 

where and (3 are constants. This Ansatz is similar to the energy density used in [2I, 
and is consistent with the EOS obtained with the approach based on the density functional 
theory j^, 4]. Let's assume that (I27p is an appropriate model for nuclear matter and that it 
satisifes the saturation condition fl22l) . This will be our model II. Once again we calculate 
the energy per nucleon E = e/ps then Taylor expand it around the equilibrium density po 
up to second order (12T|) and find 

E{pb) =ai + a2Po + ^a^po^ + a^pB^ - 2a-iPBPQ + /3VVb (28) 

Now, using (!28l) to evaluate and its derivatives we find: 

Vh = Ga^PB^pB - 4:asPo^pB + /3V(VVb) (29) 

and 

a A dpB ^ nd ,^2 N /q^^ 

In model III we consider hadronic matter at arbitrary constant baryon density, but now 
no saturation condition is imposed. This last choice is motivated by a future study of dense 
stars. In this case we calculate the enthalpy directly from fl7r|) and flTSl) obtaining the 
following expressions for the derivatives: 

Vh = Ga^PB^PB + 2^2 Vps + /3V(VVb) (31) 

and 



IV. THE KDV EQUATION 



We restrict ourselves to the one 



In this section we repeat the steps developed in 
dimensional case (a;,t) and introduce dimensionless variables for the baryon density and 
velocity: 

. PB . V 

p= — , v = — (33) 

Po Cs 
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We next define the "stretched coordinates" ^ and r as in [l|, |9|, [ill] : 



1/2 (-^ Cgt) 



r = CT 



3/2 



(34) 



where is a size scale and cr is a small (0 < a < 1) expansion parameter chosen to be 

I M - C., I 



(35) 



where u is the propagation speed of the perturbation in question. We then expand fl55]) 
around the equilibrium values: 



p = 1 + api + cr p2 + • • • 



■0 = avi + cr f 2 + 



(36) 
(37) 



After the expansion above fll6l) and fll7l) will contain power series in a (in practice we go up 
to cr^). Since the coefficients in these series are independent of each other we get a set of 
equations, which, when combined, lead to the KdV equation for pi. 



dr V2 2/iijCs2 



dpi ( ^Po \ d^pi ^ 



with the condition 



and where 



\2pB 

{<l>Po + <p)po 



1 



6^3 models II and III 

2 

^ model I 



2^2 model III 
—Aa^po model II 



-2Mcs 
PO 



model I 



u = 



13 models II and III 



St) 



model I 



(38) 



(39) 



(40) 



(41) 



(42) 



The equation fl38l) has a well known soliton solution. We may rewrite the last equation 
back in the x — t space obtaining a KdV-like equation for pi with the following analytical 
solitonic solution: 

3(u - c) 



pi(x,t) 



+ 



M 2 ^ 2mbCs2 



- c. 



sech 



l pBCs[U - Cs 

2wpo 



-{x — ut) 



(43) 



where pi = api. This solution is a bump wich propagates with speed u, without dissipation 
and preserving its shape. The expressions given by fl38|) . fl39|) and fH3|) depend on the choices 
given by (UnD, dm) and (112]). 

In model I (MQHD), the constraint (jSHD implies that Pb = M and the general equation 
becomes: 

dpi , .2^. dpi , f 9v'^Po \ ^Pi 



Illi _L _ 2n c'Pi I / Pij \ ^ Pi 

+ Cs )Pi + {2McsVRy 







with the solution given by: 



Pi{x,t) 



(3 



— sech"^ 



{u-Cs)csM 



9v 



2po 



(X 



ut) 



(44) 



(45) 



As a consitency check we take the non-relativistic limit, which, in this case, means taking a 
small sound speed 0. In this limit (3 — c^^) = 3 and (l44l) and (H5l) coincide the results 



previously obtained in 



3: 



dpi , o dpi 



2Mcs'^m^^R^) (9^3 



and 



pi{x,t) 



u 



sech 



{U - Cs)CsM 



[X 



ut) 



(46) 



(47) 



Qv V 2po 

In the limit where is large the factor 3/(3 — c^) will enhance the soliton amplitude with 
respect to the non-relativistic case. This indicates that in a medium with a stiffer EOS the 
energy propagation through solitary waves is more efficient. 

It is interesting to observe the supersonic nature of the solutions ( H5l) and ( 1471) . which is 
manifest in the arguments of the square roots. As a final remark about P3|) we notice that, 
for 



n < Cs 

2pBCs^ 2 



(48) 



the solution ( l43l) becomes negative and, in view of ( l36l) . can be interpreted as a rarefaction 
wave. A solution of this type was found in Q where nuclear matter was described by an 
EOS based on the Skyrme force. 

V. CONCLUSIONS 



The existence of KdV solitons in nuclear matter has potential applications in nuclear 
physics at intermediate energies |l| and also possibly at high energies. The experimental 
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measurements of jet quenching and related phenomena performed at RHIC [12| offer an 
unique opportunity of studying supersonic motion in hot and dense hadronic matter. With 
this scenario in mind we took the first steps in the adaptation of the KdV sohton formahsm 
to the new environment. We have extended the resuhs of our previous work [2I, showing 
that it is possible to obtain the KdV solitons in relativistic hydrodynamics. Moreover we 
have explored other equations of sate. Taking the non-relativistic limit (c^ — > 0) we were 
able to recover the previous results. 
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